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Abstract 

Using molecular dynamic simulations we study a waterlike model confined between two fixed 
hydrophobic plates. The system is tested for density, diffusion and structural anomalous behavior 
and compared with the bulk results. Within the range of confining distances we had explored 
we observe that in the pressure-temperature phase diagram the temperature of maximum density 
(TMD line), the temperature of maximum and minimum diffusion occur at lower temperatures 
when compared with the bulk values. For distances between the two layers below a certain thresh- 
old, d < dc, only two layers of particles are formed, for d > dc three or more layers are formed. 
In the case of three layers the central layer stays liquid while the contact layers crystallize. This 
result is in agreement with simulations for atomistic models. 

PACS numbers: 64.70.Pf, 82.70.Dd, 83.10.Rs, 61.20.Ja 



I. INTRODUCTION 



Water has several peculiar thermodynamic and dynamic properties not observed in other 
liquids. This is the case of the density at room pressure that has a maximum at A^C ^ ^ while 
in most materials the density increases monotonically with the decrease of the temperature. 
In addition, between 0.1 MPa and 190 MPa water also exhibits an anomalous increase of 
compressibility^ ^ and, at atmospheric pressure, an increase of isobaric heat capacity upon 
cooling^. Besides the thermodynamic anomalies water also exhibits an unusual behavior 
in its mobility. The diffusion coefficient that for normal liquids increases with the decrease 
of pressure, for water it has a maximum at 4^(7 for 1.5 atm^^. The presence of the large 
increase in the response function induced the idea of the existence of two liquid phases 
and a critical point^. This critical point is located at the supercooled region beyond the 
line of homogeneous nucleation and thus cannot be experimentally measured. In order 
to circumvent this inconvenience, experiments in confined water were performecF^. They 
showed that the large increase of the specific heat it is actually a peak that can be associated 
with the Widom line, the continuation of the coexistence line beyond a critical point. 

The drawback of experiments in nanoscale confinement is that the results obtained do 
not necessarily lead to conclusions at the bulk level. Notwithstanding this disadvantage 
the study of confined water by itself is interesting since water is present in ionic channels, 
proteins, vesicles and other cellular structures under nanoscale confinement. In order to 
understand the behavior of water under these limitations a number of experiments and 
simulations of confined water have been performed. 

Several types of confinement have been explored: experiments in cylindrical porouJ^^^^ 
and simulations in carbon nanotube j^^ * ^^ *, simulations in porous matriceJ^^'^, experi- 
ment^^nHH g^j^^j simulations in rough surfaced^^HlSI q^^i^ simulations in flat plate^^SESllESI^ 

In particular, x-ray and neutron scattering with water in nanopores show that the liquid 
state persists down to temperature much lower than in bulkP'^^'^. In these experiments in 
the case of hydrophobic walls the liquid-crystal transition occurs at lower temperatures than 
in the case of hydrophilic wall d^^ ^ ^M In some scattering experiments there are indications 
of the formation of cubic ice instead of the hexagonal ice present in the bul k^^ * ^^ i Several 
of these experiments show evidences of the presence of layers, one close to the walls and 
one at the centeij ^^ ' ^^ ' ^^ l For certain type of walls, the central crystallizes before the wall 
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layer jnini^ Therefore, the experimental results are not conclusive. They indicate that the 
crystallization in confined water depends strongly on the size of the porou^ ^^ ' ^^ ' ^^ ' ^^ ' and 
on the level of hydration water under surface^^^^. 

However, diffraction studies give only indirect information about the existence of crys- 
talline or amorphous states in water, because the Bragg peaks of ice are quite hard to distin- 
guish from liquid states. Moreover, the presence of layers is also only obtained from indirect 
evidences. In order to circumvent the difficulties of obtaining the structure of water inside 
the confined system from experiments, a number of simulations have been performe d^^ * ^^ * ^^ *. 
They employ atomistic models such as ^^id TIPSI^ and coarse-grained mod- 

gl^SESl £qj. ^ater. 

Simulations indicate that for both hydrophobi d^^ ^ ^^ ^^^^ and hydrophilicpSEZl surfaces two 
or three layers are formed depending on the distance between the confining surfaces. In 
the case of hydrophobic walls, there is a phase transition between the two to the three 
layers regime and for a certain temperature and layer separation the central layer stays 
liquid while the molecules at the walls crystallizes. In addition, in the case of hydrophobic 
walls the temperature of maximum density and the temperature of maximum and minimum 
diffusivity move to lower temperatures when compared with the bulk result^^. At very 
low pressures, cavitation appears^''. In the case of the hydrophilic walls, in agreement with 
the experiment 1 results, the system remains liquid for temperatures below the temperatures 
in the bulk cas^. 

Thermodynamic anomalies do not occur only in water, experiments for Te^^, Ga, Bi^ 
^4i|42| Gei^Te^s^^ ^ liquid metal^and graphite^^ and simulations for silicaP^, silicoii^ 
and BeF^ shown that these system also have thermodynamic anomalies. In addition, 
silicapS'^ and silicoiJ^ show diffusion anomalous behavior. In principle this systems under 
confinement could also show a shift in the anomalous properties and layering without having 
hydrogen bonds. Atomistic and coarse-grained model^SSEH for water are an interesting tool 
for understanding water and its properties, however they are not appropriated for seeking 
for universal mechanisms that would be common for water and the materials cited above in 
which the hydrogen bonds are not present but still they present the anomalous behavior of 
water. 

Acknowledging that core softened (CS) potentials may engender density and diffusion 
anomalous behavior, a number of CS potentials were proposed to model the anisotropic sys- 
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terns described above. They possess a repulsive core that exhibits a region of softening where 
the slope changes dramatically. This region can be a shoulder or a ram 

iPHSEHia. Despite 

their simplicity, such models had successfully reproduced the thermodynamic, dynamic, and 
structural anomalous behavior present in bulk liquid water. They also predict the existence 
of a second critical point hypothesized by Poole and collaborators^. This suggests that some 
of the unusual properties observed in water can be quite universal and possibly present in 
other systems. 

In this work we study the effect of the confinement in particles interacting through a 
CS potential. Our core-softened model introduced to study bulk system does not have 
any directionality and therefore it is not water. However, it does exhibit the density, the 
diffusion and the response functions anomalies observed in water. This suggests that some of 
the anomalous properties that are attributed directionality of water can be found in spherical 
symmetry systems. We explore that also some of the properties of water under confinement 
such as the presence of layering and the shift to lower temperatures of maximum density 
and of maximum and minimum of the diffusion coefficient can also be obtained with CS 
potentials. 

The paper is organized as follows: in Sec. II we introduce the model; in Sec. Ill the 
methods and simulation details are described; the results are given in Sec. IV; and finally, 
the conclusions in Sec. V. 



II. 



THE MODEL 



We study a system of particles with diameter confined between two fixed plates. 
The surfaces are formed by particles with diameter which are organized in a square lattice 
of area . The center-to-center plates distance is rf* = rf/cTp. A schematic depiction of the 
system is shown in Fig. [T] 

The particles confined between the two plates interact through an isotropic effective 
potential given by 

1 f r — Tq] 
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The first term is a standard Lennard- Jones (LJ) 12 — 6 potential with e depth plus a 
Gaussian centered on radius r = Tq and width c. We used parameters a = 5, ro/cTp = 0.7 
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FIG. 1: Model system of particles confined between plates. 



and c = 1. The pressure versus temperature phase diagram of this system in the bulk was 
studied by Oliveira et al.®'^. They found that a system or particles interacting through 
this potential exhibits a region in the pressure-temperature phase diagram where the density 
and diffusion coefficient are anomalous. 
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FIG. 2: Isotropic effective potential Eq. ([T]) of interaction between the waterlike particles. The 
potential and the distances are in dimensionless units, V — U /e and r"" — r ja^ and the parameters 
are a = 5, ro/crp = 0.7 and c—\. The inset shows a zoom in the very small atr active part of the 
potential. 



This potential has two length scales with a repulsive shoulder at r/cTp 1 and a very 
small attractive well at r jo^ ^ 3.8 (Fig. [T]). Depending of the choice of the parameters 
a, 6, c and a^, it can represent a whole family of intermolecular interactions. In this paper 
we employ a = 5, ro/cTp = 0.7 and c = 1. 
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The particle-plate interaction is given by Weeks-Chandler- Andersen Lennard- Jones po- 
tential, namelj^^'^, 



where ULj{r) is a standard 12-6 LJ potential. The cutoff distance is Vcw = S-^^^a^p, where 
(^wp = + <^i(;)/2 is the Lorentz-Berthelot mixing rule^° used when two kinds of particles 
are interacting between them. In our model, dp = (Jw = (^wp- 

III. THE METHODS AND SIMULATION DETAILS 

The system has 507 particles confined between the plates with area and distant rf, 
resulting in a number density p = N/{dL'^). The plates are located at z = and z = 
rf, whereas in x and y directions periodic boundary conditions are used. The repulsive 
interactions with the plates underestimates the number density, so we need to calculate the 
eflFective density using the effective distance de perpendicular to the plates. The new density 
will be p = N/{deL'^)^ where de ^ d — {(jp + (J^)/2 is an approach for the effective distance 
between the plated^. 

Molecular dynamics simulations at the NVT-constant ensemble and the Nose-Hoovei^^^ 
thermostat were used in order to keep fixed the temperature, with coupling parameter Q = 2. 
The interaction potential between particles has a cutoff of Tc = 3.5 and this potential was 
shifted in order to have U = at Tc- 

Several densities and temperatures are done for the following distances = d/cjp between 
the plates: 4.2, 4.8, 5.5, 6.0 and 6.3. The initial configuration of the systems were set on 
solid structure and the equilibrium states reached after 2 x 10^ steps, followed by 4 x 10^ 
simulation run. The time step was 0.002 in reduced units and the average of the physical 
quantities were get with 50 descorrelated samples. The thermodynamic stability of the 
system was checked by analyzing the dependence of parallel and perpendicular pressure on 
density namely and by the behavior of the energy after the equilibrium states. 

The thermodynamics averages in parallel and perpendicular directions to the plates are 
done employing different procedures'^. Parallel pressure, Py, is computed using the virial 
expression for the x and y directions^^ElEHll^ while the perpendicular pressure. Pi, is calcu- 
lated using two distinct methods. For systems with a strong confinement, such as rf* = 4.2 
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and 4.8, the total force perpendicular to the plates is usecP^'^, 
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For the others systems with larger distances, such as rf* = 5.5, 6.0 and 6.3, the pressure P± 
is computed through the virial expression in z directiod^. 

The dynamic of the systems was studied by lateral diffusion coefficient, L^y, related with 
the mean square displacement (MSD) from Einstein relation, 

O, ^ ..m (4) 

where ry = (x^ + y'^y^'^ is the distance between the particles parallel to the plates. 

We also studied the structure of the systems by lateral radial distribution function, g^^ (ry), 
and translational order parameter, t. We calculate the g\\{r\\) in specific regions between the 
plates, and the same for parameter t. An usual definition for is 

The 9{x) is the Heaviside function and it restricts the sum of particle pairs in the same slab 
of thickness Sz = 1. We need to compute the number of particles for each region and the 
normalization volume will be cylindrical. The g\\{T\\) is proportional to the probability of 
finding a particle at a distance r\\ from a referent particle. 
The translational order parameter is defined a d^^ ' ^^ ' ^^ ' 

t= Tl^iKO-il^e (6) 

^0 

1/2 

where ^ = r\\ps is the interparticle distance in the direction parallel to the plates scaled 
by = {Niayer/L'^y^'^. Niayer IS the avcragc of particles for each slab supposing that 
this number not change significantly (well-defined layers)^^. We use = as cutoff 

distance. 

When the system is an ideal gas, with = 1, we obtain t = 0, because the system 

is not structured. But, as the system becomes more structured, like a crystal phase, the 
5'll(r||) 7^ 1, so parameter t assumes large values. 
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All physical quantities are shown in reduced unitd^ as 
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IV. RESULTS 

Systems With Three Layers 

The first set of systems that we study corresponds to plates separated by the distances 
rf* = 5.5, 6.0 and 6.3. In aU these cases, the particles are structured in three layers in z 
direction divided in two contact layers, near to the plates, and one middle layer, located in 
the center of the plates. The formation of layering structures in confined water was also 
observed in atomistic model^^^*^. The layering density can be seen in Fig. [3]that illustrates 
the (i* = 6.3 case: (a) the snapshot of the system with T* = 0.220 and p* = 0.141, (b) the 
transversal density profile for T* = 0.220 and various densities and c) the transversal density 
profile for p* = 0.141 and different temperatures. The layers become more defined at low 
temperatures and high densities. Now we need to identify if the different layers are in the 
solid or in the liquid state. In order to answer to this question the structure is analyzed. 

The Fig. [4] shows the radial distribution function for rf* = 6.3 in two cases: (a) p* = 0.181 
and T* = 0.220 and (b) p* = 0.217 and T* = 0.140. For the case (a) the radial distribution 
of the central layer and of the contact layer are liquid-like. The contact layer shows a 
distribution compatible with a very structured liquid. For the case (b), the central layer is 
also liquid-like, however the contact layer is solid-like. The liquid-solid transition occurs at 
different temperatures and densities for the confinement rf* = 5.5, 6.0 and rf* = 6.3 that do 
exhibit three layers analyzed here. This result is in agreement with observations for SPC/E 
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FIG. 3: In (a), we have a snapshot of the system with T* = 0.220 and p* = 0.141 after the 
equihbrium state. Furthermore, the transversal density profile is shown for (b) T* = 0.220 and 
different densities, and for (c) p* = 0.141 and different temperatures. We can see the formation of 
two contact layers, near to the plates, and one middle layer. This system corresponds to d* = 6.3, 
whereas the cases like d* = 5.5 and 6.0 present the same behavior in relation to layering density. 
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FIG. 4: Radial distribution function for p* = 0.181 and T* = 0.220 in (a), and p* = 0.217 and 
T* = 0.140 in (b). Bold lines represent the ^||(^||) for the middle layer and the dashed lines 
represent the ^||(^||) for the contact layer. 

In the case of the bulkP' the potential exhibit an anomalous behavior in the translational 
order parameter t*. For normal systems the t* grows with the density, however for our CS 
potential it has a region where it does decreases with the increase of the density. Here, we 
test if this anomalous behavior is also observed in the central and contact layers. 

The Fig. [5] (a) and (b) show the translational order parameter defined by Eq.[6]as function 
of density for differents temperatures T* = 0.170,0.190,0.205,0.220,0.232,0.245,0.260 and 
0.300, at d* = 6.3, for the contact layer and for the middle layer, respectively. The dots 
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FIG. 5: Translational order parameter versus density for d* = 6.3 for (a) contact layer and the (b) 
middle layer, separately. The solid lines represent the isotherms T* = 0.170, 0.190, 0.205, 0.220, 
0.232, 0.245, 0.260 and 0.300 from the top to the bottom. The dashed lines connect the extremes. 

represent the simulation data and the solid lines identify the isotherms. The dashed lines, 
Pt-max < P* < Pt-mim identifies the region in which t* is anomalous, namely it decreases 
with the increase of the density. The region in the pressure-temperature phase diagram in 
which t* computed under confinement occurs at lower temperatures when compared with 
the bulk values^^. 

The translational order parameter for the contact layer for very low temperatures and 
very high densities is larger than the value for the middle layer what suggests that becomes 
very large the indicating the crystallization the particles at the wall are more structured 
than in the center. For the confinement with rf* = 5.5 and 6.0, a similar behavior for 

is observed. The structural differences between the layers present in our model were 
also observed in water but also in colloidal suspension, for several kinds of particle-plates 
interaction^^SMHSI confined SPC/E water by hydrophobic walpSl. 

Another propertie that exhibits anomalies in the bulk is the diffusion coefficient. In 
normal liquids, the diffusion at constant temperature grows with decreasing density, but in 
waterlike liquids there is a region < p* < Pu^^^jnax) where the diffusion decreases 

with decreasing density, so this is an anomalous behavior. In our system, this anomaly can be 
observed in the bulkP'. How the confinement affects the region in the pressure-temperature 
phase diagram where the diffusion is anomalous? In order to answer this question the lateral 
diffusion coefficient was computed as function of density as shown in Fig. |6] for (a) rf* = 5.5, 
(b) = 6.0 and (c) rf* = 6.3. In these cases, the diffusion coefficient has a region in 
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which it grows with density representing the density anomalous region. The temperature of 
maximum and minimum diffusion coefficient are lower in the confined system than in the 
bulk cas^. Our finding are in agreement with the observations of the diffusion coefficient 
in coarse-grained model confined between smooth hydrophobic plates, separated at 0.5nn]P^ 
and atomistic model^. 




p* 



FIG. 6: Lateral diffusion coefficient as function of density for (a) d* = 5.5 and isotherms T* = 0.140, 
0.160, 0.175, 0.190, 0.205, 0.220, 0.240, 0.275 and 0.320, (b) d* = 6.0 and isotherms T* = 0.140, 
0.160, 0.175, 0.190, 0.205, 0.220, 0.235, 0.250 and 0.290 and (c) d* = 6.3 and isotherms T* = 0.140, 
0.170, 0.190, 0.205, 0.232, 0.260 and 0.300, from the bottom to the top. The dots represent the 
simulation data and the solid line is just a polinomial fit to isotherms. The dashed lines connect 
the extremes of diffusion. 

In the bulk, our potential exhibits a density anomalous region in the pressure-temperature 
phase diagram. How confining affects the TMD line? In order to answer to this question 
the TMD is computed under confinement. The density anomaly is given by {dp/dT)p^^ = 0, 
so using the following Maxwell relation 

(i),,=-(^),-0/ 

density anomaly can be found through {dP\\/dT)p = 0, what is equivalent to the minimum 
of the parallel pressure versus temperature. 

Fig. [T] illustrates the parallel pressure versus temperature phase diagram for the (a) 
d* = 5.5, (b) d* = 6.0 and (c) d* = 6.3 cases. The thin solid lines are the isochores, the 
solid bold lines represent the TMD, the dashed lines are the lateral diffusion extremes, the 
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dashed-dotted lines are the translational order parameter extremes for a contact layer and 
the dotted lines are the translational order parameter extremes for the middle layer. 
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FIG. 7: Phase diagrams P* -T* for (a) d* = 5.5, (b) d* = 6.0 and (c) d* = 6.3. The thin solid lines 
are the isochores (a) 0.117 < p* < 0.233, (b) 0.105 < p* < 0.209 and (c) 0.099 < p* < 0.253. The 
solid bold line represents the TMD, the dashed line are the lateral diffusion extremes, the dashed- 
dotted line are the translational order parameter extremes for a contact layer and the dotted line 
are the translational order parameter extremes for the middle layer. 

A comparison between the TMD of the confined and the bulk systems is given by Fig. |8] 
(a). The densities and pressures ranges of the TMD's locations are also shown in Table [ij 
The TMD lines of confined systems are shifted to lower temperatures and higher densities 
when compared with the bulk. 

For d* = 5.5 and 6.0, the TMD lines are shifted to higher pressures, whereas for d* = 6.3 
this shifting occurs to slightly lower pressures when compared with the bulk TMD. The 
non monotonic shift in pressure when compared with the bulk results can be attibuted to 
the fact that we employ the lateral pressure for the confined system while we use the total 
pressure for the bulk system. 

Kumar et al.^^ found that the TMD line for confined systems is shifted to lower tempera- 
tures but in the same range of pressures when compared with the bulk system . For TIP4P 
water model in contact with six hydrophobic spheres, Gallo and Rovere^^ found that TMD 
line in confined systems is shifted to lower temperatures and higher pressures. Furthermore, 
they observed that the spinodal curve follow the shifting of the TMD line. Similar result is 
observed in our systems. Xu and MolinercP^, using a coarse-grained model for water (mW, 
Monatomic Water Model)^^^ confined in nanopores of diameter 1.5nm, also found a TMD 
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FIG. 8: Phase diagram P|* — T* comparing in (a) the shifting of the TMD hnes for the confined 
systems in relation to the bulk. A comparison between the anomalies of the confined systems is 
given in (b) for diffusion and in (c) for translational order parameter of the middle layer. For (b) 
and (c), the solid bold lines represent the system with d* = 5.5, the dashed-dotted lines represent 
(i* = 6.0 and the thin dotted lines represent = 6.3. 

shifted to lower temperatures and higher pressures. So, the density anomalies observed in 
our systems have a good agreement with other atomistic and coarse-grained simulations. 

TABLE I: Densities and pressures ranges of TMD's location of confined and bulk systems. 



d* 


density range 


pressure range 


5.5 


0.149 <p* < 0.188 


0.800 < < 1.237 


6.0 


0.129 <p* < 0.162 


0.599 < < 0.965 


6.3 


0.122 < p* < 0.153 


0.541 < P* < 0.875 


bulk 


0.110 < p* < 0.140 


0.552 <P* < 0.913 



The perpendicular pressure is shown in Fig. [9] as funtion of (a) temperature and (b) 
density, at rf* = 6.3. A monotonic increasing behavior is observed in both cases. The 
ranges of densities and temperatures are the same done in P|* — T* phase diagram. The 
other systems (rf* = 5.5 and 6.0) have a similar behavior and they are not shown here 
for simplicity. The same behavior for perpendicular pressure also was observed by Kumar 
et al.ES. 
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FIG. 9: Phase diagrams (a) - T* and (b) - p* at d* = 6.3. The dots are the simulation 
data and the sohd hne just connect the isochores in (a) and the isotherms in (b). The ranges of 
densities and temperatures are the same done in P|* — T* phase diagram. 

System With Two Layers 



The confinement by very narrow distances induces the transition from three to two layers. 
In this subsection, we study a system with plates separated by rf* = 4.2. A snapshot in Fig.[TOl 



(a) shows the two contact layers, without middle layer. Figs. 10 (b) illustrate the behavior 
of the transversal density profile for fixed temperature, T* = 0.220, but for diflFerent total 
densities. Figs. 10 (b) also shows the density profile but for fixed total density, p* = 0.155, 



and several temperatures. The struct ur at ion in just two contact layers is due the strong effect 
of confinement. Structure of bilayer is observed for hydrophobic confinement in TIP5P'^2''^, 
TIP4F'^ and mW^ models of water. 
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FIG. 10: (a) Snapshot of the system after the equilibrium state. In (b) we can see the tranver- 
sal density profile for T* = 0.220 and several densities, and in (c) for p* = 0.155 and several 
temperatures. 
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Fig. 11 (a) shows the parallel pressure versus temperature phase diagram. The isochores 
are represented by thin lines, and the TMD by the solid bold line. It also shows the diffusion 
extrema (dashed line) and the translational order parameter extrema (dashed-dotted line). 
The comparison between the TMD line for the rf* = 4.2 case and the TMD for the bulk 
system is illustrated in Fig. |lT| (b). The location of the TMD is 0.604 < P|[ < 0.959 and 
0.137 < p* < 0.170. Similarly to what happens in the three layer system, the TMD line 
shifts to lower temperatures when compared with the bulk. This characteristic is evidenced 
in Fig.[lT](c). 
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FIG. 11: Phase diagrams P|* — T* showing the density, dynamic and structural anomalies in (a). 
In (b) we have a comparison between the TMD lines of d* = 4.2 and the bulk system, and in (c) 
the same comparison is done with TMD lines for systems that present formation of three layers. 



Fig. 
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(a) illustrates the perpendicular pressure versus temperature phase diagram for 
various isochores that show no TMD line. No anomalous behavior is observed similarly to 
what happens in the three layers regime. Fig. [12] also shows the radial distribution function 
(in (b)) and the mean square displacement (in (c)) for the lateral direction. For p* = 0.137, 
these figures show a amorphous solid-like behavior for T* = 0.160 and a liquid-like behavior 
for T* = 0.270. A solid-to-liquid transition was also observed for the TIP5P modeW^^ for 
the TIP4P modeJ^ and for the mW model^. 

An anomalous region for translational order parameter and for lateral diflFusion are also 
observed in this extremely confined system as can be seen in Figs. [13] (a) and (b), respectively. 
The dashed lines connect the extremes of these anomalies, defining the anomalous region 
that we can see in P|* — T* phase diagram in Fig. jlTj (a). 
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FIG. 12: In (a), the monotonic increasing behavior is observed for perpendicular pressure with 
the temperature. In (b), we have a ^||(^||) for a amorphous sohdhke state, with T* = 0.160, and a 
hquidhke state, with T* = 0.270, both at p* = 0.137. The MSD is observed for these cases in (c). 
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FIG. 13: In (a), the translational order parameter as function of density, and in (b) the lateral 
diffusion coefficient as function of density. The isotherms are connecting the simulation data 
(points) and they are from the top to the bottom T* = 0.175, 0.190, 0.205, 0.220, 0.235, 0.250, 
0.270 and 0.320 for t* (a) and from the bottom to the top for D*^ (b). The dashed lines connect 
the extremes observed in both cases. 

Three-to-Two layers System 



The system with rf* = 4.8 exhibits an unusual behavior that resembles properties of both 



two layers and the three layers systems. In Fig. 14 (a) the lateral diffusion as function 



of density is illustrated. The diffusion anomaly is only present at very low temperatures, 



T* < 0.150. In addition, as shown in the Fig. 14 (b) the TMD line is also restricted to low 



temperatures. In Fig. 14 (c) the comparison between the TMD line of the confined system 
with (i* = 4.8 and the TMD of bulk confirms this shift of the TMD to very low T*, much 
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lower than the shift observed for confinement with 4.8 < rf* and 4.8 > rf* discussed above. 
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FIG. 14: In (a), the lateral diffusion coefficient as function of density for isotherms from the 
bottom to the top T* = 0.118, 0.130, 0.140, 0.150, 0.160, 0.190 and 0.220. Diffusion extremes 
are connected by dashed lines. In (b), we have the phase diagram showing the TMD line and the 
diffusion extremes. A comparison between the TMD line of this system with bulk is given in (c). 



This (i* = 4.8 case is not only peculiar for exhibiting a very low temperature of maximum 
density but also for presenting anomalous behavior at the perpendicular pressure versus 



temperature phase diagram. In Fig. 15 (a) the isochores at the perpendicular pressure 
versus temperature phase diagram exhibits minima not shown for 4.8 > rf* > 4.8. In order 
to shade some light in the reason for the unusual behavior of the system for rf* = 4.8 we 
explore the behavior of the structure and of the stability of the layers. Fig. [15] (b) shows 
the lateral pressure versus density of fixed temperatures. For T* < 0.190, the slope of the 
curve first increases and then decreases. This change even thougth not indicating a phase 
transition is usually observed before the phase separation would be stablishecP^. 

The Fig. [Ts] (c) shows the tranversal density profile for p* = 0.139 and for a number of 
temperatures. At this density, three layers are present for low temperatures, T* = 0.118 and 
0.150. At high temperatures, T* = 0.250, two layers are well defined with particles equally 
distributed between them without forming a third layer. 

In order to understand the two-to-three layers transition, the change in the structure 
of the central layer with temperature and density is checked. The structure of the middle 



layer is shown in the Fig. [16] for different temperatures and densities. At low temperatures 
(cases (a) at T* = 0.118 and (b) at T* = 0.160) and high densities, p* > 0.170, the central 
layer is divided in many sublayers. By decreasing the density the many layers give rise at 
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T 



P 



FIG. 15: The perpendicular pressure versus temperature is given in (a) and versus density is 
given in (b). The simulation data (dots) are connected by solid lines for better visualization. The 
transversal density profile for p* = 0.139 and some temperatures is shown in (c). 

p* = 0.170 to a central layer that disappears as the density is decreased any further. At 
high temperatures, T* > 0.160 (case (c) shown at T* = 0.220), as the density decreases the 
system passes from three-to-two layers without forming the sublayers. At low temperatures, 
due to the presence of the many sublayers the density anomalous region appears. The TMD 
originates from particles moving from one scale to the othei'^. 




z z z 

FIG. 16: Tranversal density profile for the middle layer at many densities and (a) T* = 0.118, (b) 
T* = 0.160 and (c) T* = 0.220. The dashed lines represent states for < 0.170, the solid lines 
for p* > 0.170 and the bold circles are for p* = 0.170. 

Our interpretation is in accordance with the assumption of Kumar et al.^, that the split 
of the middle layer in sublayers is justified by the density anomaly. Using SPC/E model 
for water and hydrophobic rough plates, Giovambattista et alW^ found a phase transition 
between a bilayer ice and a trilayer heterogeneous fluid for different distances between the 
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plates. For smooth plates separated at 0.8nm and SPC/E model, Lombardo et al.'^ also 
found a phase transition between two and three layers. 

Changes in diffusion coefficient as function of separation between the plates are reported 
in systems with transitions between 6 and 5, 5 and 4, 4 and 3 layer J^^, suggesting that the 
dynamic behavior change in systems with structural transitions. So, the dynamic behavior 
of our system at rf* = 4.8 is another possible explanation for the peculiar behavior on its 
anomalies. 

V. CONCLUSIONS 

In this paper, the effects of confining a system of particles interacting through core- 
softened is explored. 

The formation of three layers, two close to the walls and one central is observed for 
large values of rf*, while two layers are observed for small values of rf*. In addition the 
region in the pressure-temperature phase diagram where the density anomaly appears moves 
to lower temperatures. These results are similar to the results obtained in atomistic and 
coarse-grained models where unlike our model the directionality of the h-bonds is explicitly 
includecp™. 

Our results indicate that layers are formed in order to minimize the particle-particle 
interaction potential and the wall-particle interaction. Therefore, if the walls are distance 

= 5.0 it is possible to fit three layers in which each will is distant 2.0 from the other 
and the contact layer is distant 1.0 from the wall. For rf* = 4.0 it is only possible to fit 
two layers distant 2.0 from each other and 1.0 from the wall. The density, diffusion and 
structural anomalous behavior that implies particles moving from one length scales to the 
other (moving from the length scale at 3.0 to the length scales at 1.0) occurs only along 
the parallel plane, therefore the anomalies appear as function of P\\ . 

The case rf* = 4.8 is the boundary between the two layer and the three layer cases. This 
case allow us to observe how the presence of the density anomaly is related with moving 
from different particle-particle distances. 

Our results suggest that effective spherical symmetric two length scales potentials are an 
interesting tool for understanding the mechanisms that arise from confining systems with 
density, diffusion and structural anomalies. Due to their simplicity the results obtained can 
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be generalized to other experimental realizations besides water. 
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